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We present a detailed de Haas van Alphen effect study of the perovskite CaVOs, offering an un- 
precedented test of electronic structure calculations in a 3d transition metal oxide. Our experimental 
and calculated Fermi surfaces are in good agreement - but only if we ignore large orthorhombic dis- 
tortions of the cubic perovskite structure. Subtle discrepancies may shed light on an apparent 
conflict between the low energy properties of CaVOs, which are those of a simple metal, and high 
energy probes which reveal strong correlations that place CaVOs on the verge of a metal-insulator 
transition. 

PACS numbers: 71.18.+y, 71.20.-b, 71.30.+h 



Transition metal oxides (TMOs) have for many 
decades been a rich source of novel and intriguing elec- 
tronic behaviour. The narrow bands in the TMOs are 
particularly susceptible to correlation effects, leading to 
phenomena as unexpected and diverse as high tempera- 
ture superconductivity [|l[ |[, colossal magnetoresistance 
[3|, and the recently discovered "orbitons" in LaMnOa 
[4, ||. The Mott transition - a metal-insulator 
transition (MIT) driven by electron-electron interactions 
- lies at the heart of TMO physics, and the 3(i^ sys- 
tem CaVOs is emerging as an important model com- 
pound that lies just on the metallic side of this transition 

While the general principles involved in the MIT are 
probably understood, the predictive power of existing 
theories is poor when applied to real materials. On 
the one hand, models of the MIT that explicitly tackle 
dynamical many-body interactions use generic Hamilto- 
nians that ignore the effects of real crystal structures; 
on the other, single-particle self-consistent band struc- 
ture calculations based on the Local Density Approxi- 
mation (LDA), which have had the most success at pre- 
dicting the electronic structure of real metals, fail spec- 
tacularly in insulating TMOs such as NiO which is pre- 
dicted by the LDA to be metallic. This problem is com- 
pounded by a lack of high quality data on the electronic 
structure of TMOs: for example until now the exclu- 
sive source of information on the Fermi surface geometry 
in TMOs has been angle resolved photoemission spec- 
troscopy (ARPES), but this technique detects a much 
higher energy response than is observed in transport and 
thermodynamic measurements, and the results are fre- 
quently distorted by surface effects. 

CaVOs is one of a series of 3d^ perovskites that strad- 
dle the MIT. It is a metal, but large orthorhombic dis- 
tortions of the cubic perovskite crystal structure and the 
spectroscopic properties of CaVOs suggest that it lies 
very close to the MIT [|[ ^, 10, ll|]. Photoemission spec- 
troscopy in CaVOs appears to show: (a) the formation of 



Hubbard bands [|| |T|, |T|, || — a precursor of the insu- 
lating state arising from strong electron correlation; and 
(b) a weaker spectral intensity at the Fermi energy Ey 
than is predicted by band structure calculations [|l2[ ^ . 
From these and other measurements a broad picture is 
emerging of a striking conflict between low and high en- 
ergy scales in CaVOs, the former being revealed for ex- 
ample in specific heat measurements and suggesting that 
CaVOs is a simple metal with small mass enhancements 
(normally a signature of weak electron-electron interac- 
tions), while the latter are probed by spectroscopic meth- 
ods and show strong electron-electron interactions. The 
origin of this conflict is unclear. It may be an experi- 
mental artifact arising from enhanced sensitivity of spec- 
troscopy to the surface as opposed to the bulk, but extra 
care has been taken to minimise this An intriguing 
suggestion is that the conflict, and the breakdown of the 
LDA that it implies, can be explained by a strong inter- 
action that is non-local in space but local in time, such as 
for example the unscreened Coulomb interaction. But for 
this to occur, we require either a breakdown of screening 
in a material that at dc-frequencies at least has a normal 
metallic conductivity, or else a new non-local interaction 
such as might be mediated by orbital fluctuations. 

Further experimental data are clearly needed, and this 
was the motivation for the present de Haas van Alphen 
(dHvA) effect study, which allows an unprecedented and 
detailed comparison between modern electronic structure 
calculations and an experimentally determined Fermi 
surface in a 3d TMO. It should be noted that LDA cal- 
culations have an excellent record predicting the Fermi 
surface shapes of even strongly correlated electron sys- 
tems such as UPts 0, in which electron-electron in- 
teractions are enormously strong. This is probably be- 
cause the Fermi volume is invariant under interactions 
]T8[ : volume conserving distortions of the Fermi surface 
require not only a /c-dependent (i.e. spatially non-local) 
interaction such as was invoked to explain the pho- 
toemission results, but also an anisotropic Fermi surface. 
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FIG. 1: Example of a dHvA field sweep on CaVOs with the 
applied field parallel to the 6-axis. The lower inset is a blow-up 
of the high-field region. The upper inset shows a logarithmic 
plot of the dHvA frequency spectrum, taken over the 18 T to 
15 T field range. 



so that states in different regions of the Fermi surface 
undergo different shifts in energy. 

Quantum oscillations, such as the magnetisation oscil- 
lations in the de Haas- van Alphen (dHvA) effect, arise 
from the quantisation of cyclotron motion of charge car- 
riers in a magnetic field. This provides the most detailed 
and the most reliable information about Fermi surface 
properties [|o|. A prerequisite for the observation of the 
quantum oscillation effects is that quasiparticles must be 
able to complete a cyclotron orbit without scattering — 
this typically requires mean-free paths of ~ 1000 A or 
more. It is this requirement of very high purity single 
crystals that has hitherto prevented quantum oscillations 
from being observed in 3d TMOs. 

We used high-quality single crystals of CaVOa, grown 
by the floating zone method yielding single-crystalline 
grains of typical size ~ 4 x 4 x 6 mm^ . We kept the oxy- 
gen off-stoichiometry as small as possible by annealing 
the sample in air at ~ 150°C for four days [^. The off- 
stoichiometry of Ca and V was below the detectability 
limit (less than 1 %) of an inductively coupled plasma 
atomic emission spectrometer. The CaVOa crystals are 
orthorhombic (Pnma) containing four CaVOa formulae 
in the unit cell with lattice parameters a = 5.314 A, 
b = 7.521 A, c = 5.339 A (Nakotte, H. & Jung, M.-H., 
private communication); the samples were oriented by 
Laue x-ray diffraction. 

We have performed a thorough dHvA rotation study on 
two high-quality single crystals of CaVOa (1x1x2 mm^, 
residual resistivity po c:^ 1.5 /il^cm) with the direction of 
the applied magnetic field B rotating in the {a-b) plane 
from b to a, and in the (a-c) plane from a to c. The 
experiments were carried out in a low-noise supercon- 
ducting magnet system in field sweeps from 18 T to 10 T, 



at temperatures down to 50 mK. A modulation field of 
14.3 mT amplitude was applied to the sample, and the 
second harmonic of the voltage induced at the pick-up 
coil around the sample was recorded, essentially measur- 
ing d'^M/dB'^. 

A typical signal trace, demonstrating the high qual- 
ity of our data, can be viewed in Fig. ^ The Fourier 
transform as a function of 1/B gives the dHvA fre- 
quency spectrum, with peaks at a frequency F corre- 
sponding to extremal (i.e., maximum or minimum) FS 
cross-sectional areas S measured in planes perpendicular 
to B, via S = 27ieF/n. 

In Fig. ^ we present a density plot of the dHvA spec- 
trum as a function of the direction of 5, representing the 
information of about 100 field sweeps at angular intervals 
of 2°. In comparison, we have calculated the dHvA fre- 
quencies Fcaic from four FS sheets (see Fig. ^) obtained 
by a band calculation using the computer code KANSAI- 
94 with the LDA scheme for the one-electron 

exchange-correlation potential. In the calculation, we 
used basis functions in |^ + G| < i^max = 6.20 (27r/A), 
where k is the wavevector and G is a reciprocal lattice 
vector, to obtain about 1440 basis linear augmented plane 
waves; the core and valence states were calculated self- 
consistently by the scalar-relativistic scheme with 
muffin-tin radii of 0.25 a for Ca, 0.20 a for V, and 0.14 a 
for O. 

As seen in Fig. |^, there are very significant discrepan- 
cies between Fcaic (thin dashed lines) and the experimen- 
tal dHvA frequencies Fexp (grey-shaded). To list a few, 
(a) the ~ 9kT branches in Fexp have no correspondence 
in Fcaic, and (b) the ^63 sheet is not observed in Fexp un- 
less we accept a deviation of about 20%. Only ^64 and 
some of the orbits on ^62 show satisfactory coincidence. 

Although we observe this large discrepancy between 
experiment and LDA, we do not believe that this is the 
end of the story, because we have found that excellent 
agreement between theory and experiment can be re- 
gained if we ignore the orthorhombic distortion in the 
crystal structure^ and recalculate the Fermi surface as 
though CaVOa were cubic. We attempted to fit the cal- 
culated LDA bands of the orthorhombic lattice structure 
to a hypothetical "cubic" tight-binding (TB) Hamilto- 
nian consisting of V 3dxy^ 3dyz and 3dzx orbitals, with 
nearest neighbour, next-nearest neighbour and hybridiza- 
tion matrix elements. Then, the results (thick solid lines 
in Fig. H) show extremely good agreement with exper- 
iment: the differences in F are only 4% or less. As 
depicted in Fig. ^, in this assumption, the orthorhom- 
bic FS sheets ^61, ^62, and ^63 will coalesce into the 
hypothetical "cubic" FS sheets P and 7, while the ^61 
sheet remains unaffected. The 7 sheet connects up with 
itself in the neighbouring BZ, forming a "jungle gym"- 
type structure. The large hole orbit in this "jungle gym" 
explains the high-frequency features in Fexp, e.g. along 
the 6-axis. 
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FIG. 2: Density plot of the dHvA spectrum of CaVOs as a function of the direction of B (greyshaded) . The main panel uses 
the high-field (18 T 15 T) sections of each sweep, while the lower inset uses all available data (18 T 10 T), with higher 
frequency resolution at the expense of signal-to-noise. The right inset shows a blow-up around the a-axis for fields in the 
ac-plane. Superposed are the frequency predictions from LDA band structure calculations, for both the orthorhombic (thin 
dashed lines) and the hypothetical cubic (thick solid lines) unit cells — the colour code corresponds to Fig. ^ 



We also deduced the effective quasiparticle masses m* 
from the thermal broadening factor Rt = X/ sinhX, 
X = 27r'^ kBTm* /eh by comparing dHvA amplitudes 
at different temperatures in runs from 50 mK up to 1.1 K. 
As seen in Fig. ^ we have m* ex F for all orbits and for 
all field directions assessed, save for the hole orbit of the 
"jungle gym". Such homogeneous scaling of m* with F 
is seen in the TB model in two-dimensional metals with 
a low-filling and with primarily nearest-neighbour inter- 
actions. This results also supports our assumption of 
the effective cubic symmetry, because the FS sheets of a 
cubic system are essentially formed from three inter- 
secting cylinders containing the dxy^ dxz^ and dyz states 
respectively, and each has a two-dimensional character. 
The specific heat deduced from dHvA is 8.6mJ/moleK, 
which is in good agreement with the bulk experimental 
value of 7.3mJ/moleK and is enhanced over the un- 
renormalized LDA density of states by about a factor of 
two. 

Although some details in the experimental dHvA spec- 
trum of CaVOa remain yet to be identified and are possi- 
bly associated with remnant features of the intrinsic or- 
thorhombic structure, the overall angle dependence of the 



dHvA frequencies is altogether well captured by the hy- 
pothetical "cubic" TB model, rather than the orthorhom- 
bic LDA band calculation. 

Our apparently ad hoc assumption that the effective 
crystal symmetry sensed by the ?>d electrons is cubic, 
may be justified if the gaps in the energy bands that 
are produced by the orthorhombic distortion are small 
enough — an effect known as magnetic breakdown (MB) 
[|o|. However, its appearance here seems to be surpris- 
ing, since the orthorhombic distortion in CaVOa is not 
small (the buckling of the V-O-V bond angle is 160°). 

We hope that our results will inspire further, more de- 
tailed electronic structure calculations. It will be inter- 
esting to see if a high-resolution LDA calculation will 
show the small gaps due to the orthorhombic distortion 
that our results imply, and thereby reproduce the predic- 
tions of our tight binding model. If not, some more exotic 
explanation of the non-appearance of orthorhombic gaps 
may be required. Also, while the differences between our 
tight-binding Fermi surface and the measurements are 
small, they are systematic. In particular, it can be seen 
in the lower insert of Fig. ^ that the fit to the (3 surface is 
clearly worse than the others. The (3 sheet is composed 
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the FS sheets and thus the electronic structures near the 
Mott transition in CaVOa. The experimental FS cannot 
be explained by the LDA band calculation for the intrin- 
sic orthorhombic lattice, but the results are well captured 
by a hypothetical "cubic" electronic structure. At low en- 
ergy scales, CaVOa is seen to be a standard metal with 
moderate mass enhancement, in marked contrast to the 
strong correlations visible at higher energies. While fur- 
ther efforts are required to comprehend this discrepancy 
and to gain a quantitative understanding of the origin of 
the MB effects, this work provides at least a first glimpse 
of the completely unexplored fermiology near the Mott 
transition. 



FIG. 3: FS sheets predicted by LDA. The right panels show 
the FS sheets in a cubic TB fit to the LDA data, while on the 
left side these sheets are folded over into the orthorhombic 
BZ. The field directions in the rotation study are indicated 
around the surfaces in the top row. 
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FIG. 4: The ratio of the effective masses to the free electron 
mass for all observed Fexp in CaVOs at various field orien- 
tations, spaced by about 15° along b ^ a ^ c. The masses 
scale with the dHvA frequency, except for the 9kT hole orbit. 



of states that are all near an intersection of the three 
cylinders of the cubic Fermi surface model, so perhaps 
this is indeed a signature of a non-local interactions, but 
detailed computational work is clearly required to clarify 
this, starting with a high-resolution LDA calculation of 
the electronic structure. 

In summary, this first comprehensive investigation of 



We thank D. M. Broun, H. Bando, and W. Y. Liang 
for their support. We also thank H. Nakotte and M.-H. 
Jung for sharing their data prior to publication. This 
work was supported in part by the JST Overseas Re- 
search Fellowship and the UK EPSRC. 



[1] J. G. Bednorz and K. A. Miiller, Z. Phys. B 64, 189 
(1986). 

[2] J. Orenstein and A. J. Millis, Science 288, 468 (2000). 

[3] Y. Tokura (ed.). Colossal Magnetoresistive Oxides, Gor- 
don and Breach, London (2000). 

[4] E. Saitoh et a/.. Nature 410, 180 (2000). 

[5] Y. Tokura and N. Nagaosa, Science 288, 462 (2000). 

[6] N. F. Mott, Metal- Insulator Transitions, Taylor and 
Francis, London (1990). 

[7] M. Imada, A. Fujimori and Y. Tokura, Rev. Mod. Phys. 
70, 1039 (1998). 

[8] I. H. Inoue et a/., Physica B 194 1067 (1994). 

[9] A. Fukushima et a/., J. Phys. Soc. Jpn. 63 409 (1994). 
[10] N. Shirakawa et aL, J. Phys. Soc. Jpn. 64 4824 (1995). 
[11] I. H. Inoue et a/., Phys. Rev. B 58, 4372 (1998). 
[12] I. H. Inoue et a/., Phys. Rev. Lett. 74, 2539 (1995). 
[13] K. Morikawa et al, Phys. Rev. B 52, 13 711 (1995). 
[14] M. J. Rozenberg et a/., Phys. Rev. Lett. 76, 4781 (1996). 
[15] H. Makino et a/., Phys. Rev. B 58, 4384 (1998). 
[16] K. Maiti et al, Europhys. Lett. 55, 246 (2001). 
[17] L. Taillefer and G. G. Lonzarich, Phys. Rev. Lett. 60, 
1570 (1988). 

[18] J. M. Luttinger, Phys. Rev. 119, 1153 (1960). 

[19] More technically, the self-energy ^{k,co) has to be k- 
dependent. This is usually ignored in theories of the MIT, 
but it seems plausible that as screening breaks down on 
the verge of the MIT such a /c-dependence could arise. 

[20] D. Shoenberg, Magnetic Oscillations in Metals, Cam- 
bridge University Press, Cambridge (1984). 

[21] O. Gunnarsson and B. I. Lundqvist, Phys. Rev. B 13, 
4274 (1976). 

[22] D. D. Koelling and B. N. Harmon, J. Phys. C: Solid State 

Phys. 10, 3107 (1977). 
[23] I. H. Inoue, Researches of the Electrotechni- 
cal Laboratory 986, 1 (1999). (available at 



http:// www.etl.go.ip/ en / results / researches / researches986- 



the qHvA effect in a 3d TMO has unraveled the shape of 



^.html[ 



